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ABSTRACT 

We consider a distributed sensor network in which the sin- 
gle antenna sensor nodes observe a deterministic unknown 
parameter and after encoding the observed signal with a phase 
parameter, the sensor nodes transmit it simultaneously to a 
multi-antenna fusion center (FC). The FC optimizes the phase 
encoding parameter and feeds it back to the sensor nodes such 
that the variance of estimation error can be minimized. We 
relax the phase optimization problem to a semidefmite pro- 
gramming problem and the numerical results show that the 
performance of the proposed method is close to the theoretical 
bound. Also, asymptotic results show that when the number 
of sensors is very large and the variance of the distance be- 
tween the sensor nodes and FC is small, multiple antennas do 
not provide a benefit compared with a single antenna system; 
when the number of antennas M is large and the measurement 
noise at the sensor nodes is small compared with the additive 
noise at the FC, the estimation error variance can be reduced 
by a factor of M. 

Index Terms — Distributed sensor network, multi-antenna 
fusion center, maximum likelihood estimation, phase-only 
analog encoding, asymptotic estimation error. 

1. INTRODUCTION 

Sensor networks have been widely studied for detection and 
estimation problems. Recently, considerable research has 
focused on the fusion of analog data in a distributed sen- 
sor network to improve estimation performance. In (TJ, the 
authors considered estimation problems involving single an- 
tenna sensors and a single antenna FC. The sensor nodes use 
an amplify-and-forward scheme to transmit observations to 
the FC over fading channels and the optimal allocation of 
power to the sensors was investigated. In [2], the asymptotic 
variance of the best linear unbiased estimator for a single 
antenna distributed sensor network was derived and the effect 
of phase quantization errors were analyzed. Four different 
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multiple access schemes for a single antenna decentralized 
estimation system with a Gaussian source were investigated 
in Q, and the scaling laws for large number of sensors were 
derived. A coherent multiple access channel was considered 
in 10 and the optimal linear decentralized estimation scheme 
was investigated. It is well known that for a multiple an- 
tenna communications system, the link capacity generally 
increases linearly with the minimum number of antennas at 
the transmitter or receiver. It is also expected that the estima- 
tion performance of a sensor network would also benefit from 
a multi-antenna FC, although prior work on this scenario is 
limited. A system with a multi-antenna FC was considered 
in 0, which showed that for Rayleigh fading channels, the 
reduction in estimation error variance is bounded by 2 when 
the number of sensors approaches infinity. 

In this paper we consider a distributed sensor network in 
which several single antenna sensor nodes observe a deter- 
ministic parameter corrupted by noise, and simultaneously 
transmit the observed signal to a multi-antenna FC. The low- 
complexity sensor nodes are assumed to transmit an analog 
signal with constant power and adjustable phase. The FC de- 
termines the optimal value of the phase for each sensor in 
order to minimize the maximum likelihood (ML) estimation 
error, and the FC feeds this information back to the sensors so 
that they can encode their observed signal accordingly. The 
estimation performance of the phase-optimized sensor net- 
work is shown to be considerably improved compared with 
sensors that use non-optimized phase. 

The paper is organized as follows. Section [2] describes 
the system model and Section [3] formulates the phase opti- 
mization problem and proposes a numerical solution based 
on semidefinite programming (SDP). Asymptotic analyses for 
large number of sensors and large number of antennas are pro- 
vided in Section [4] Numerical results are then presented in 
Section|5]and our conclusions are summarized in Section|6] 

2. SYSTEM MODEL 

We assume that N single-antenna sensors in a distributed sen- 
sor network independently observe a deterministic parameter 
8 6 C. The sensor nodes encode the observed signal with a 
phase parameter a; (1 < i < N) and simultaneously transmit 



it to the FC. Assuming the FC is configured with M antennas, 
the received signal at the FC is expressed as 

y = Ha6» + HDv + n, 

where H = [hi , . . . , hjv] and h^ e C A/X 1 is the channel vec- 
tor between the zth sensor and the FC, a = {ai, . . . , ajv} T 
contains the adjustable phase parameters (|a.;| = 1), D = 
diagjai, . . . , a^r}, v is the measurement noise at the sen- 
sor nodes which is assumed to be Gaussian distributed with 
covariance V = E{vv H } = diag {0^,1 > • • ■ ,o^jv}> n 
is the additive Gaussian noise at the FC with covariance 
E{nn H } = crfJ-M and Im is an M x M identity matrix. 

Assuming the FC is aware of the channel matrix H, the 
noise covariance V and cr 2 , the FC calculates the ML estimate 
of using |6] 



'ML — 



a g H g (HVH fl + all M )- l y 
a H H H (HVH ff + cr2i M )-iHa' 



The estimator Oml is unbiased and the variance of Oml is 
given by 

Var(lW) = (a H H H (HVH H + ^i M )-i H a) _1 . (1) 
The variance is lower bounded by 

Var(0 Mi ) > M — (HH(HV ^ - a 2 lM) - lu) . (2) 
where A max (-) denotes the largest eigenvalue of a matrix. 

3. PHASE-ONLY ANALOG ENCODING METHOD 

From (0, we see that Var(6>|H)ML is inversely proportional 
to a quadratic form in a, which suggests the possibility of 
minimizing the variance by adjusting the phase of the sig- 
nals transmitted by the sensors. The advantage of adjusting 
only the phase is that such adjustments do not impact the co- 
variance of the sensor noise observed at the FC, unlike what 
would occur if the transmit power of the sensors was also ad- 
justed. This is the reason for the assumption of sensors with 
constant transmit power. 

We formulate the following optimization problem at the 

FC: 



min Var(#Mi) 

a 

s.t. \cii\ = 1, i = 1, . . . , N. 



(3) 



Define B = H fl (HVH ff + o£Im) -1 H, then we can rewrite 
problem (0 as 



max a^Ba 

a 

s.t. \aA = 1, i = 1, . . . , N. 



(4) 



If there are only two sensors in the network, a closed- 
form solution to problem (0 can be obtained. Defining 



B 



c be jl3 
be-^ c 
value of B is given by 



with b, c > 0, then the largest eigen- 



Amax(B) = C+b, 

and the solution to (0]i is given by the corresponding eigen- 
vector 



a = \e 3h 



JI32 



where fli — 02 — ft- 

For the general situation where N > 2, we relax prob- 
lem (O and convert it to a standard SDP problem that can be 
solved with standard tools such as cvx [7|. To begin with, 
we rewrite (O into the following equivalent form 



max 

A 

S.t. 



tr(BA) (5) 

A M = 1, i = l,...,N 
rank(A) = 1 
A y 0, 



where Aj^ denotes (i, i)th element of A. Relaxing the rank- 
one constraint on A, we can convert problem (0 to a standard 
SDP problem: 



max 

A 

S.t. 



tr(BA) 
A >- 0. 



(6) 



.N 



Defining B r = real{B}, B^ = imag{B}, and similarly for 
A,, and Aj, © can be converted to the equivalent real form 



max tr(B r A,. — B^A; 

{A r ,Aj} 



(7) 



s.t. 



A r j 5 j — 1, i — 1 
A i A. r 



.N 



>- 0. 



Problem (0 can be efficiently solved using a standard interior 
point method. Denote the solution to problem (0 as A*. If 
rank(A*) > 1, we can use a method similar to Algorithm 2 
in |0] to extract a rank-one solution. In Section the numer- 
ical results show that the performance of the proposed phase 
encoding method is close to the theoretical lower bound in (0. 

4. ASYMPTOTIC ANALYSIS 

In the asymptotic analyses conducted in this section, we as- 
sume a path-loss based model 

hi = — hj, 

df 

where di denotes the distance between the ith sensor and the 
FC, a is the path loss exponent and h; denotes the normal- 
ized channel component. Furthermore, the elements of the 



normalized channel are assumed to have constant amplitude 
and random phase: 



hi = [e j "t^,e> 



. ,e 



3H,m-\T 



where 7.^ is assumed to be uniformly distributed over [0, 2n]. 

4.1. Performance Bound for Large Number of Sensors 

From (O, the lower bound of V&y(0\H)ml depends on the 
largest eigenvalue of U H (HVH H + oIImY 1 **- To eval- 
uate this eigenvalue, we approximate HVH H as a diagonal 
matrix. The (m, n)th element of HVH fl is given by 



HVH* „ = £ 



W V.l 



According to the strong law of large numbers, when N — > oo, 
the following equation holds 



1 N 

i=l 



(b) I E 

I 0, 



^}e{^.-^.")} 



m 7^ n, 



where (a) follows from the assumption that "fi >m , d{ and a 2 i 
are independent of each other and (b) is due to the fact that 
7i im and 7i in are independent and uniformly distributed over 
[0, 2ir}. Thus, when N — > oo, we can approximate HVH fl 
as 



lim HVH H 

N^>oo 



NE • 



where E |^S^-| = J2iLi Based on (O, we have 

-A max (H H H) 



A max (H ff (HVH fl + ^I M )"iH) 
1 



NE 



(10) 



where (c) is due to A max (H ff H) = A max (HH ff ). Plugging 
([Tol l into (fj), the following lower bound for Nsi(9ml) can be 
obtained: 



lim Var( 

N^oo 



'ML, 



> 



NT N - 



(11) 



Clearly, an upper bound on the estimate variance can be 
found by considering the single-antenna case, which results 
in 

(12) 



lim Var(6»ML) < 

N—too 



(t n ± 



When N — > oo, the ratio of the lower to the upper bound is 
given by 

(gkiT (■{*»' . Mil 



lim - „ ' - = - — ) ' J ' =1 ' J . (13) 

— e{^} e{^} 



Interestingly, we see from ( fT3l that when Var { ^sr j <C 

E | |, the gap between the bounds is small and the avail- 
ability of multiple antennas at the FC does not provide a 
benefit compared with the single antenna system. 

4.2. Scaling Law with Number of Antennas 

Using the matrix inversion lemma, the inverse of noise co- 
variance matrix can be written as 

H H (HVH H + all M )- 1 H 



H 



H 



* H ( V-V^H^H 



n H I h 



(14) 



The (m, n)th element of H H H is given by 



M 



1 

u m u n i=1 



J'(7r»,i-7m,i) 



(9) Similar to {§), when M — > oo, the following holds: 



M 

lim -V 

M^oo M ^ 



I . m = n, 
0, m ^ n. 



Based on (Q2), we can approximate H H H as 



lim H ff H«Mdiag<f-J-,- 

Substituting ( TTol ) into ( TBI , we have 
H H (HVH H + ct 2 Im )-i H 

M 



' r/ 2Q 



(15) 



(16) 



diag ■ 



M 



d{ a a 2 n + Ma 



2 ' 
v.i 



Since H if (HVH fl + o- 2 Im) _1 H is thus approximately di- 
agonal, for any phase vector a, we have 



M 



lim Var (9 M l) 



(17) 



From (TPD ). it can be observed that when Ma% 4 <C cr^df , a 
reduction in the estimate variance by a factor of M can ap- 
proximately be achieved, and to realize this gain, the phase 
vector a can be selected arbitrarily. 
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Fig. 1. Performance of the ML estimator with increasing 
number of sensors for M = 4. 

5. NUMERICAL RESULTS 

To evaluate the performance of the proposed approach, sev- 
eral numerical experiments were carried out. In the numerical 
examples, the transmit power of sensor nodes are normalized 
to 1 and the path loss exponent a is set to 1. Both dj and 
al i are assumed to be uniformly distributed, with di ranging 
from 2 to 7, i between 0.001 to 0.01 and the additive noise 
power <r„ is set to 0.1. In the results, each point on the curve is 
obtained by averaging over 300 realizations of h^. To evaluate 
the performance without feedback, a is set to a vector of all 
ones (the random phase component is subsumed in the chan- 
nel vector). The results in Fig. 1 show that the performance 
of the proposed method is very close to the lower bound in 
eq. (f2]) and significantly better than without feedback. Fig. 2 
shows that the benefit of the feedback reduces as the number 
of FC antennas grows, although the gain is still significant for 
reasonable values of M below 10. We see in both figures that 
as either N or M get large, our asymptotic expressions match 
the numerical results well. 

6. CONCLUSION 

In this paper, we investigated a phase-only analog encoding 
scheme for a distributed sensor network composed of single- 
antenna sensors and a multi-antenna FC. We relaxed the phase 
optimization problem to an SDP and the numerical results 
show that the performance of the proposed method is close 
to the theoretical lower bound. Also, asymptotic results for 
cases with large numbers of sensors or antennas were de- 
rived. The asymptotic results indicate that when the number 
of senors is large and the variance of the distance between the 
sensor nodes and the FC is small, the availability of multiple 
antennas at the FC does not provide a benefit compared with 
the single antenna case. On the other hand, when the number 
of antennas M is large and the measurement noise and the 



Fig. 2. Performance of the ML estimator with increasing 
number of antennas for N = 4. 

additive noise satisfy Ma% t <C a^d\ a , a reduction in estima- 
tion variance of M can be obtained. 
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